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^D ■ The DC Josephson effect is theoretically studied in a planar junction in 

which a two dimensional electron gas (2DEG) infinite in lateral directions is 
■^ I in contact with two superconducting electrodes placed on top of the 2DEG. 

An energy gap in the excitation spectrum is created in the 2DEG due to 
C^ ' the proximity effect. It is shown that under certain conditions, the region of 

the 2DEG underneath the superconductors is analogous to a superconducting 
— ^ , region with an order parameter egexp(i^), where Eg {eg < A) depends on 

Ch ', the interface transmittance and the Fermi velocity mismatch between the 

TT \ superconductors and the 2DEG. 
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I. INTRODUCTION 

In recent years significant progress has been made in the preparation and study of Joseph- 
son junctions, in which the weak couphng is reahzed through a two dimensional electron gas 
(2DEG). These systems are prepared using semiconductor heterostructures with supercon- 
ducting contacts. The structure Nb-InAs(2DEG)-Nb appears to be the most promisingpQ, 
although it may be possible to obtain the Josephson effect also in other structures based, 
for example, on GaAsoEi. In Josephson junctions with a 2DEG as a weak link, one expects 
phenomena analogous to the conductance quantization in quantum point contactsETa. One 
of these phenomena is the predicted quantization of the critical Josephson current /JlSlliJ. In 
short quantum contacts [L < C,, C, = hvp/irA is the coherence length), the critical current 
is predicted to increase with increasing contact width in a step like way, with a step height 
equal to eA/h. 

The geometry of the contacts shown in Fig. la has been considered in previous theo- 
retical workslldO. In this case, electrons move in a channel or 2D region of finite length 
and experience Andreev reflections at the SN interfaces. The excitation spectrum of the 
2DEG (or IDEG) is changed due to the interference of waves reflecting at the opposite SN 
boundaries; in particular, bound states decaying into the S regions and corresponding to 
energies e < A appear in the systemll3TJ. These bound states give the main contribution 
to the Josephson current Ic at low temperatures. 

It is, because of its practical significance, of interest to study the DC Josephson effect 
in the geometry shown in Fig. lb. In this geometry, the 2DEG is unbounded in the {x, y) 
plane (the z axis is directed perpendicular to the 2DEG plane). We will find the condensate 
Green's functions in the N layer, i.e., in the layer with the 2DEG, and show that due to 
the proximity effect, the properties of the N layer at |a;| > L/2, where L is the distance 
between both superconductors, are analogous to those of a superconductor with an effective 
order parameter egexp{i(j)) with Eg dependent on the interface transmittance and the Fermi 
velocity mismatch. We note that the pair potential in the N layer An determined from the 
selfconsistency equation equals zero because the electron-phonon coupling constant in the N 
layer is supposed to be vanishing. For a general discussion of the proximity effect in terms 



of a tunneling model see Ref. |T^. We adopt the following simplified model of the system 
shown in Fig. lb. The transmittance of the SN interface Tb[x) is proposed to be dependent 
on the coordinate x; it varies from a given value Tb at |x| > L/2 to zero at |x| < L/2 on a 
characteristic length xq. We assume that xq is much larger than the Fermi wave length kp^ 
in the N layer, but is smaller than the coherence length in the 2DEG .^20 = hvp/TcSg. Then, 
we can use an adiabatic approximation to calculate the energy spectrum of the 2DEG at 
|x| > L/2. 

II. THE PROXIMITY EFFECT IN THE 2DEG. 

Let us consider the energy-diagram for a system shown schematically in Fig. 2. The 
region at 2; < is occupied by a superconductor, and there is a quantum well with the 
2DEG in the layer < z < d. The Fermi momenta in the S and N regions (pp = Jpfi + p^ 



and kp = Jkl + ky, respectively) differ greatly from each other due to a significant difference 



in electron concentrations in these regions; namely, pp = {Sn'^nsY^^ ^ fcp = (2v™2d)^^^ (the 
formula for kp is written for the case when only the lowest subband in the quantum well 
is filled). Generally speaking, there exists a Schottky barrier at 2; = 0, which we model by 
a potential of the form Ubab6{z), where Ub and a^ are the height and width of the barrier 
varying smoothly at |x| = L/2. 

In order to calculate the excitation spectrum and wave functions of the 2DEG, we write 
the well known Bogoliubov-de Gennes equations 



H^k = £fc^fc. 



(1) 



Here ^^ is the two component wave function, Sk is the excitation energy, and H is the 
Hamiltonian of the system, 
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Here mi and m,2 are the effective masses in the S and N regions respectively, and f/gm is the 
difference between the potentials of the conduction band edges in the S and N regions (see 
Fig. 2). 

One can show that propagating (into the S region) states with e^ > A and bound states 
with Ek < ^ exist. The latter states correspond to a branch of the spectrum with an effective 
energy gap Eg, which is small if the barrier transmittance is small. For low temperatures 
(T <C A), the main contribution to the critical current Ic originates from these bound modes. 
We will restrict ourselves to the case of low temperatures. A solution describing the bound 
states has the from 



\^fc (r) = exp(iA;||r|| 
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Consider first the wave function in the quantum well (the second term in Eq.(^). The 
first two terms {E±) describe an electron excitation moving forward and backward and the 
second term {H±) correspond to a hole excitation. Momenta of these excitations along the z 
axis can be presented in the form k^ = kn + Skn and kz = kn + Skn, where fc„ = {'7i/d){n + 1) 
and \6kn,Skn\ -C /c„. Here we will consider the case when only the lowest subband, with an 
energy Eq = 7r/(i^?7i2, is filled {n = 0) and other subbands (n > 1) are empty. So we are 
interested in the states with n = 0. A relation between 6kQ and 6kQ can be obtained from 
Eq.(g), 



£k = ^k + VoSko = -(6 + VoSko 



(4) 



Here C,k = {kj — kpQ)/2m,2 is the kinetic energy of electrons in the 2DEG relative to the 
Fermi energy, Ep = (/cq + fc|Q)/2m2 + f^sm? ^fo is the Fermi momentum in the limit of infinite 
barrier height, Vq = ko/m2. 



Consider now the wave functions in the S region, which decay over the length k ^. The 
coefficients {u, v) have the usual form, 

«' = ^(1 + W^p), v^ = \{l-ije,), (5) 

but in this case the functions E,p are purely imaginary 

ip = -i^Pz/mi, Ep = (A^ + g)'/'. (6) 

The momentum p^ approximately coincides with the Fermi momentum pp where we assume 
the reflection at the interface to be specular, i.e. p\\ = k\\ and, as noted above, k\\ <C pp. 
Hence for the Fermi energy, we can write 

ep = pl/2mi = (pI + pj)/2mi ^ pl/2m,. (7) 

In order to find the excitation spectrum, i.e. the dependence of e^ on the momenta in the 
{x,y) plane k\\, we must use boundary conditions at the interface {z = 0). These conditions 
consist in continuity of "^^ (z) at z = and in a relationship between the derivatives dz^^k^^ (z) 
at ^ = 00 

^kA+0) = ^kA-O), -^dAkAz)\+o - -^dz^kAz)\-o = Ubab^kM- (8) 

Substituting Eq.(^) into Eq.(|]), we get a set of algebraic equations for the coefficients E±, 
H± and B±. The solvability condition results in a dispersion relation for the excitation 
spectrum at e^ < A 

1 + akOtkiw'^ + s^) + w{ak + ak) + is{ek/ip){ak - a^) = (9) 

Here ak = Skod -C 1 and ak = Skod = —{ak + 2C,k/£o)- ui = {2Ubab'm2/ko) is a dimensionless 
parameter characterizing the barrier transmittance, the factor s = ppm2/kQmi depends 
on the mismatch of the Fermi momenta and the effective masses. Eq.(P) determines the 
spectrum of bound states with energies Sk < A. Eq.(P) can be rewritten in the form 



1 + ^''' 



(A2- 4)1/2 



W .o _2 



{^k--egor + e' (10) 



s 



where {C,k — —^go) is the relative kinetic energy of electrons moving in the {x,y) plane. The 
quantity Ego = Eq s/iw"^ + s^) is the energy gap in the excitation spectrum of the 2DEG 
induced by the proximity effect in the case of very low barrier transmittances. Indeed, under 
the condition 

(^o/A)^^^ « 1 (11) 

one can neglect the second term in the square brackets in Eq. ([To|) and obtain for Ek not too 
close to A 
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Therefore, the dependence of Ek on k\\ is nearly the same as in a 2D superconductor with 
the energy gap 

where Eq is the subband energy for n = 0. This dependence is shown in Fig. 3 for several 
values of w. In Fig. 3 it is shown that, with changing w, not only the value of Eg is changed, 
but also the k\\ at which the minimum in Ek occurs. This may be understood intuitively, 
since increasing the SN-barrier transparency will alter the exact form of the electron-hole 
wave function in the 2DEG, "^kui^)- Since kp is increased as compared to fcpo when the 
barrier transmittance is finite, the minimum in e^ is expected to occur at larger fcy. 

With increasing the temperature T the energy gap in the S region A(T) is diminished. 



and the condition Eq.(|Tl]) is violated at T sufficiently close to Tc. The dependence of the 
energy gap in the 2DEG Eg on A(T) is determined from Eq.(|IO|) if we put (^^ — 7^30) = 0, 
which means putting Ek at a minimum. Then we obtain the equation for Eg 
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This dependence is shown in Fig. 4. The maximal value of A(T)/Ego equals A{0)/Ego. If this 
value is very large, the energy gap in the excitation spectrum of the 2DEG coincides with 
Ego in a wide range of T. The characteristic temperature T* determining a transition of Eg 
from Ego to A(T) is given by the equation A(r*) = EgQ. At low temperatures (A(T) = A(0)) 
we can also calculate the influence of the barrier, expressed in w, from Eq. ([13|) . Fig. 5 shows 
the energy gap Eg in the 2DEG as function of the transmittance T^ = l/((w/2s)^ + 1) of the 
S-2DEG interface, for different values of s. 

If the excitation energy Ek exceeds A, the wave functions in the S region do not decay, 
but oscillate, and the k vector runs over a continuous set of values. These wave functions 
describe the propagation of two electrons (incident on the barrier and reflected from it) and 
a hole that appears as a result of Andreev reflection. The bound states obtained above are 
closely related to those studied earlier in the 3D caselllrEj. 

The wave functions corresponding to the bound states are determined by Eq.(^. By 
introducing new variables 

w 

£k = £k/£0 = Oik + ik/^G, tk = ik/£Q — t-w, t-w = ~^^, 2' (^^) 



tf^ + S^ ' 



we can write for the coefficients from Eq.(|^) the relations 



E- = E+exp{2ia), H^ = H^exp{2ia) 

B+ = 2ik [{E+Up + H+Vp)Ek - {E+Up - VpH+){tk + t^)] /{6^ - Ekf'' 

5_ = 2Ek \{E^Vp + R^UpYk - (E+Vp - UpH+){tk + t^)] /(<5' - Ek)'^^. (15) 

Using these relations and Eq.(|^) for the wave functions, we can find the Green's functions 
qR{A) £qj, ^YiQ 2DEG in the system shown in Fig. lb and calculate the critical Josephson 
current. 



III. THE GREEN'S FUNCTIONS AND THE DC JOSEPHSON EFFECT 



Consider the system shown in Fig. lb. Electrons move in the quantum well in the 
{x, y) plane. As shown before, the wave functions in the 2DEG change drastically due to 
the proximity effect. In particular, the condensate Green's functions F^^^^ are induced in 
the quantum well, and therefore the Josephson effect is possible in this system. In order to 
determine the critical current^ we need to know the functions F^^^\ that is, the nondiagonal 
elements of the matrix G^^ : F^'^'^^z, z';q) = G12 {z,z';q). We are interested in the 
current / averaged over the quantum well width, which means that we must find F^'^'^\k) = 
{F^^^\z, z'] q)), where the brackets denote averaging over z {0 < z < d). As is well known, 
the matrix components of G^^^^ are expressed through the components of the wave functions 
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where the sum is taken over the two branches of the spectrum determined by Eq.(|T( 
We suppose that the transmittance of the SN barrier is small, implying e^ ^ A, 
conditions s/w <^ 1 and s/w"^ <^ A/^o are fulfilled. Then, one can show that E^ 
are coupled by the relation 



and the 
and H. 



(17) 



Taking into account Eqs.(|T^ and (jT^), one can find from Eq.([TB|) that the components of 
the Green's functions are equal to 
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Here Ek 
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+ eltlf' 



is the excitation energy, e 



g is the energy gap in the excitation 
These functions are identical to the 



spectrum in the 2DEG at |a:;| > L/2 (see Eq.(|T^)) 
Green's functions of an ordinary two dimensional superconductor with a spatially dependent 
energy gap because, as supposed, the parameter w{x) varies from a constant value w at |a;| > 
L/2 to 00 at \x\ < L/2. The characteristic length of the w{x) variation is xq, which is small 
as compared to the coherence length in the 2DEG, ^20 = hvp/iiEg and large as compared 
to the Fermi wave length k^^ = (m2fF)^^- Therefore the system under consideration is 
equivalent to a 2DEG contacting at |x| > L/2 the superconducting 2DEG with the effective 
order parameter Egexp{±i(f)/2), where e^ -C A (strictly speaking this magnitude is achieved 
at |a;| > L/2 over distances of the order of the coherence length as it takes place in ordinary 
SNS junctions). Hence for the critical current Ic one can use the formulae obtained in 
Refs. 0-|l4| where the width of the N layer was assumed to be large or in Ref. 10 where a 
one dimensional channel is analyzed. If the width, Ly, of a channel in the 2DEG shown in 
Fig. lb is comparable with the Fermi wave length kp^ at |x| < L/2 and is much larger at 
|x| > L/2, we can use the expression for Ic obtained in Ref. ^ 



N{eEg/h) sin(0/2) tanh((£g/2T) cos(0/2)). 



(19) 



where N is the number of subbands below the Fermi leveL This means that Ic will increase 
step wise with increasing the electron density in the quantum well, and the height of the 
steps equal eSg/h. 

IV. CONCLUSIONS 

We have found the critical Josephson current Ic for a 2DEG in contact with two su- 
perconductors. In contrast to previously analyzed systems the effective order parameter in 
the 2DEG egexp{i(f)) is reduced in comparison with Aexp(i0) in the superconductor, and 
its magnitude is determined by the SN barrier (Schottky barrier) transmittance and by a 
mismatch of the Fermi momenta and the effective masses in the S and N regions (here 
is the macroscopic phase of the superconductors). The barrier transmittance characterized 
by the parameter w may depend on the carrier density in the 2DEG, n2D- Therefore, the 
critical current Ic will depend on n2D even when only the lowest subband is filled. 
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FIGURES 



FIG. 1. Schematic representation of S-2DEG-S Josephson junctions, as studied in Refs. 1C,11 
(la) and in the present paper (lb). The spacing between the superconductors is L, the width of 
the 2DEG is d. 

FIG. 2. Energy diagram of the system under consideration - a superconductor and a 2DEG are 
separated by a potential barrier Uf,. Hatched areas denote states filled with electrons (we suppose 
that only the lowest subband in the quantum well, with energy £q, is occupied). The energy gap 
in the excitation spectrum of the 2DEG is induced due to the presence of the superconductor. 

FIG. 3. The normalized excitation energy in the 2DEG, £k/^, vs the longitudinal momentum 
of electrons, fcy, at different values of the parameter Sgo/A = (eo/A)s/(s^ + t(j^): Ego/A = 0.2 
(dashed line); 0.3 (solid line); 0.6 (dotted line). 

FIG. 4. The normalized energy gap in the excitation spectrum of the 2DEG £g/£go vs the 
normalized energy gap in the S region A/e^o- When the temperature is increased towards Tc, 
A(T) goes to 0, thus violating the condition of Eq.([Tl|). 



FIG. 5. The normalized energy gap e^/A vs the transparency Th, expressed in the dimensionless 
parameter w, of the barrier between superconductor and 2DEG, for different values of the mismatch, 
s = p-pm2 / kQ-nii: s = 0.5 (solid line); 1 (dashed line); 2 (dotted line). 



